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The uniform electron gas (UEG) at finite temperature is of key relevance for many applications in 
dense plasmas, warm dense matter, laser excited solids and much more. Accurate thermodynamic 
data for the UEG are an essential ingredient for many-body theories, in particular, density functional 
theory. Recently, first-principle restricted path integral Monte Carlo results became available which, 
however, due to the fermion sign problem, had to be restricted to moderate degeneracy, i.e. low to 
moderate densities with = f/as -J 1. Here we present novel first-principle configuration PIMG 
results for electrons for < 1. We also present quantum statistical data within the e"*-approximation 
that are in good agreement with the simulations at small to moderate r^. 

PACS numbers: 05.30-d, 05.30.Fk, 71.10.Ca 


Thermodynamic properties of quantum degenerate 
electrons are vital for the description of matter at high 
densities, e.g. CHS], such as dense plasmas in compact 
stars or planet cores, e.g. HH, as well as in laser fu¬ 
sion experiments at NIF, e.g. Rochester HO] or 

Sandia Besides, the electron component is of 

crucial importance for understanding the properties of 
atoms, molecules and existing and novel materials. The 
most successful approach has been density functional 
theory (DFT)-combined with an approximation for the 
exchange-correlation potential. Its success is based on 
the availability of accurate zero temperature data for the 
UEG which is obtained from analytically known limit¬ 
ing cases combined with first-principle quantum Monte 
Garlo data [T5] . 

In recent years more and more applications have 
emerged where the electrons are highly excited, e.g. by 
compression of the material or by electromagnetic radi¬ 
ation (see above), which require to go beyond zero tem¬ 
perature DFT. This has led to an urgent need for accu¬ 
rate thermodynamic data of the UEG at finite temper¬ 
ature. One known limiting case is the highly degener¬ 
ate ideal Fermi gas (IFG), and perturbation theory re¬ 
sults around the IFG, starting with the Hartree-Fock 
and first order correlation corrections (Montroll-Ward) 
[H ng, are long known; for an analytical fit at high 
densities, see Ref. m and for further improved approxi¬ 
mations, such as the e^-approximation, see Refs. [IlHIS]. 
These approximations break down when the Coulomb 
interaction energy among the electrons becomes compa¬ 
rable to their kinetic energy, requiring computer simu¬ 
lations such as path integral Monte Carlo (PIMG), e.g. 
|20| . While restricted PIMG (RPIMC) results for dense 
multi-component quantum plasmas, e.g. mull] as well 
as direct fermionic PIMG (DPIMC) results have 

been available already for 15 years, only recently finite 
temperature RPIMC results for the UEG (or jellium) 
have been obtained m- It is well known that fermionic 
PIMG simulation in continuous space suffer from the 
fermion sign problem (FSP) which is known to be NP 
hard [28j . This means, with increasing quantum degen¬ 


eracy, i.e. increasing parameter y = which is the 

product of density and thermal DeBroglie wave length, 
= h?[2TTmkBT\~^, the simulations suffer an expo¬ 
nential loss of accuracy. RPIMC formally avoids the FSP 
by an additional assumption on the nodes of the den¬ 
sity matrix, however, it also cannot access high densi¬ 
ties |29j . rg < 1 [r^ = f/as, where r is the mean in¬ 
terparticle distance, = 47rf^/3 and as the Bohr ra¬ 
dius]. Also, the quality of the simulations around rg = 1, 
at low temperatures 0 = ksT/Ep < 0.125 [Ep is the 
Fermi energy] is unknown. However, this leaves out the 
high-density range that is of high importance, e.g. for 
deuterium-tritium implosions at NIF where mass densi¬ 
ties of 400 gcm“^ (up to 1596 gcm“^) have recently been 
reported [5] (are expected along the implosion path 0), 
corresponding to « 0.24 (rg = 0.15), see Fig. 

The authors of Ref. ^7} also performed DPIMC sim¬ 
ulations which confirmed that, for 0 < 0.5 and rg < 4, 
these simulations are practically not possible, see Fig. 
We also mention independent recent DPIMC simula¬ 
tions [30] that are overall in good agreement with the 
data of Ref. m but indicate large deviations for the low¬ 
est temperatures and rg ;< 2 (lower energies). Finally, 
a recent attempt to avoid the FSP by an approximate 
treatment of exchange cycles m also yielded lower en¬ 
ergies for Tg = 1 and 0 = 0.125 than the RPIMC data of 
Ref. [27]. To bridge the gap between the known analyti¬ 
cal result for the ideal Fermi gas at rg = 0 and previous 
simulations (rg > 1) and to provide comprehensive in¬ 
put data for finite temperature DFT, several fits have 
been proposed recently [33] [33] . However, they crucially 
depend on the quality of the underlying simulation data. 

It is the purpose of this Letter to improve this situ¬ 
ation. We present the first ab initio simulation results 
that avoid a simplified treatment of fermionic exchange 
for rg < 1 and finite temperatures, 0 < 1.0. We apply 
the recently developed fermionic configuration path in¬ 
tegral Monte Carlo (CPIMC) approach to the UEG and 
demonstrate its capabilities for 33 spin polarized elec¬ 
trons in a cubic box of side length L (as was studied 
in [3I1I33]). Our simulations have no sign problem for 
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density parameter r. 


FIG. 1: Density-temperature plain in the warm dense matter 
range. IGF: typical inertial confinement fusion parameters 
[ 8 ]. Quantum (classical) behavior dominates below (above) 
the line 0 = 1 . F = e^/ffcsT is the classical coupling pa¬ 
rameter. Red dots: available finite temperature RPIMG m 
and DPIMG |30| data for the UEG. Blue dots: ground state 
data of Ref. |38| . Green dots: GPIMG and analytical results 
of this work. 


0 < Ts < 0.4 and are accurate up to = 1. 

CPIMC for the UEG. The jellium Hamiltonian in 
second quantization with respect to plane waves {r\k) = 
_^^ik r ^ w G has the familiar form 

(we use Rydberg units throughout) 

^ = ^ ^ + Em, (1) 

i i<j,k<l 

— dvre ^ki-\-kj ^kk-\-ki 

"^ijkl ~ '^ijkl ~ Wijlk, Wijkl — ^ , ( 2 ) 

^ik 

where the first (second) term describes the kinetic (inter¬ 
action) energy, and kik = ki — kk- The Madelung energy 
Em accounts for the self-interaction of the Ewald sum¬ 
mation in periodic boundary conditions [34] for which we 
found Em ~ —2.837297 • {?,/. The operator 
a\ (Sj) creates (annihilates) a particle in the orbital \ki). 
In the interaction term, the ki = k^ and kj = ki com¬ 
ponents cancel with the interactions with the positive 
background. While the complete (infinite) set of plane 
waves (r \ki) forms a basis in the single-particle Hilbert 
space, for simulations it has to be truncated at a number 
Nb of orbitals. 

In conventional RPIMC and DPIMC, the system Q 
is treated in the coordinate representation allowing for 
a numerically exact description in the classical strongly 
coupled limit and for weak degeneracy. CPIMC [3S], in 
contrast, is constructed in a way that it allows for ex¬ 
act simulations in the opposite limit of the ideal Fermi 
gas, Tg = 0 [36], and at weak to moderate coupling 
and strong degeneracy. This is achieved by represent¬ 
ing the iV-electron state in second quantization m as 


Si = (2, 5, 0,3) |{n(2)}) = 1001110...) 


5 

4 

13 

° 1 
0 


Tl T2 T3 

imaginary time r 


FIG. 2: Typical closed path in Slater determinant (Fock) 
space. The state with three occupied orbitals |feo), |fei), Ifea) 
undergoes a two-particle excitation (si,ri) that replaces the 
occupied orbitals \ko),\k 3 ) by |fc 2 ),|A: 5 ). Two further excita¬ 
tions occur at T 2 and T 3 . The states at the “imaginary times” 
T = 0 and T = P coincide. All possible paths contribute to 
the partition function Z, Eq. ©• 


|{n}) = |ni, 712,...), where the are fermionic occupa¬ 
tion numbers (rii = 0,1) of the orbitals \ki). This way, 
fermionic anti-symmetry is “built in” exactly. The parti¬ 
tion function Z and quantum-statistical expectation val¬ 
ues, such as the internal energy U, are straightforwardly 
computed in Fock space as 

Z(0,rg;IV)=Tr|{„pe-'3^, (3) 

U{e,rs-,N) = {H) = Z-iTr|{„piJe-'5", (4) 

avoiding the sum over A^! permutations over the canon¬ 
ical density operator that is the origin of the no¬ 

torious FSP in DPIMC. The trace is evaluated using the 
concept of continuous time PIMC which has been suc¬ 
cessfully applied to bosonic lattice models We 

have generalized this concept to continuous fermionic sys¬ 
tems with long range interactions [33 HI] . The main idea 
is to split the Hamiltonian into a diagonal, D, and an off- 
diagonal part, Y and summing up the entire perturbation 
series of the density operator in terms of Y. The 

final result, for the case of the UEG, is [IS]: 
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O. = 

Si ^pqrs 
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l—p l—r 


(5) 


( 6 ) 


where Si = {p,q,r,s), with p < q, r < s denotes a 
quadruple of pairwise different orbital indices. 

Thus the partition function is represented as a sum 
over /3-periodic “paths” in Fock space, in imaginary time. 
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which we illustrate in Fig. For an ideal Fermi sys¬ 
tem a path is characterized by a single 7V-particle Slater 
determinant |{n}). For a correlated Fermi system the 
original determinant |{n}) = (straight horizon¬ 

tal lines in Fig.[^ is interrupted by excitations of the type 
(s,r): at time r, a pair of occupied orbitals \kr)i |A:s) is 
replaced by the previously empty pair \kp), \kq). Paths 
differ by the number K of excitations (“kinks”), their 
times Ti... Tfe on the r-interval [0, /3] and the involved 
quadruples of orbitals Si... sk- The partition function 
clearly reflects this summation over the different types 
of kinks, integration over the kink times and summation 
over K [cf. first line of Eq. ([^] . The weight of each path 
[terms in the second line of Eq. ([^] is determined by the 
Fock state matrix elements of the Hamiltonian, where di¬ 
agonal elements Eq. arise from the kinetic 

energy and the mean-field parr of the Coulomb interac¬ 
tion, whereas off-diagonal elements, (—are due 
to the remaining Coulomb interaction (correlation part) 
pS] . Expression ([^ is exact for Nb —>■ oo allowing for ab 
initio thermodynamic simulations of the UEG. 

Formula (§ and similar expressions for thermody¬ 
namic observables, such as the internal energy [45) . can 
be efficiently evaluated using Metropolis Monte Carlo. 
To this end, we developed an algorithm that generates 
all possible paths in Slater determinant space thereby 
assuring ergodicity. For the UEG a total of 6 different 
steps are required, including addition and removal of a 
single kink and pairs of kinks, modification of an exist¬ 
ing kink and excitation of single orbitals. The complete 
set of steps and the associated MC algorithm cannot be 
presented here, for details see Ref. [IS]. 

Numerical results. Our finite-temperature CPIMC 
algorithm was extensively tested for Coulomb interact¬ 
ing fermions in a ID harmonic trap [35] . A first test of 
the present algorithm for the UEG for iV = 4 particles 
showed excellent agreement with exact diagonalization 
data [IS], and the results were confirmed by density ma¬ 
trix QMC IH]. Here, we extend these simulations to 
N — 33 particles. First we check the convergence with 
respect to the basis size Nb and show a typical case in 
Fig. [3|a, for = 0.4 and four temperatures. The seal- 
ing with respect to a: = 1 /Ng ' |47| allows for a reliable 
extrapolation to a; —)■ 0 and to set Nb to 2109 for all sim¬ 
ulations, giving a relative basis incompleteness error not 
exceeding the statistical error (Icr standard deviation). 

With these parameters, we have performed extensive 
CPIMC simulations in a broad range of r^- and 0-values. 
In contrast to DPIMC and RPIMC, our ab initio simu¬ 
lations (without any simplifications besides the choice of 
Nb) pose no problem for the ideal and weakly coupled 
UEG, up to Tg ~ 0.4. For larger r^, we observe a rapid 
decrease of the average sign, in analogy to the harmonic 
oscillator case [35]. This gives rise to convergence prob¬ 
lems of the MC algorithm in case a path with many kinks 
is attempted. WeMherefore, introduce an artificial kink 
potential in Eq. |^, I4(iF) = -p 1]“^, for 
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FIG. 3: Convergence of the CPIMC simulations, a. Con¬ 
vergence with the single-particle basis size Nb for Vs = 0.4. 
The expected scaling with m is well reproduced, b. 

Convergence with respect to the kink potential parameter k 
(see text) and extrapolation to 1/k —> 0, corresponding to 
K —>■ oo, for Vs = 1.0 and 6 — 0.0625. The asymptotic value 
is enclosed between the red and blue line. 


calculations with rg > 0.4, yielding the correct parti¬ 
tion function in the limit k —> oo. Performing simula¬ 
tions for different a, we generally observe a rapid conver¬ 
gence of the total energy allowing for an extrapolation 
to !/«;—>■ 0. This is demonstrated for the most difficult 
case (rg = 1,0 = 0.0625) in Fig. §b. The asymptotic 
value and the error estimate are computed from the two 
extreme cases of a horizontal and linear extrapolation. 
With this procedure the simulations could be extended 
to Tg = 1, with the total error not exceeding 0.15%. 

Our simulation results for the exchange-correlation en¬ 
ergy per electron E^c are summarized in Fig.[^ The data 
cover the whole range 0 < Cg < 1 and include the ideal 
Fermi gas where f const (Hartree-Fock limit). 

A detailed table of the various energy contributions is 
presented in the supplementary material [50]. A non¬ 
trivial observation is the non-monotonic temperature de¬ 
pendence (cf. crossing of the red and pink curves) that is 
in agreement RPIMC and the macroscopic fit of Ref. [32] . 
Interestingly, all curves seem to cross over smoothly into 
the RPIMC data [571 1121 , for Vg > 4, as indicated by 
the dotted line. There is an obvious mismatch in the 
range rg ~ I... 4. Since our curves are accurate within 
the given error, this discrepancy is expected to be due 
to the (unknown) systematic error involved in RPIMC. 
Also, the energy obtained within the novel approach of 
DuBois et al. ED, cf. the data point for rg = 1 and 
0 = 0.125 in Fig. ^ is found to be too low. 

Macroscopic results. Predictions for a macroscopic 
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FIG. 4: (Color online) Exchange-correlation energy (times r^) 
for 33 spin-polarized electrons and four temperatures. Com¬ 
parison of our CPIMC results (full symbols with error bars 
|51|1 and RPIMC results of Ref. [23 (open symbols). The dot¬ 
ted line is an interpolation between the CPIMC and RPIMC 
data for 0 = 0.0625. Also shown is the data point of DuBois 
et al. 131] for 0 = 0.125. 


system, based on data for just 33 particles, will inevitably 
lead to a loss of accuracy. Brown et al. have attempted 
such a mapping of their RPIMC data for iV = 33 to 
IV —oo and published finite size corrections (FSC) in 
the supplement of Ref. m- A ground state formula [FSC 
(a)] has been presented by Drummond et al. [52]. We 
tested both FSC, after incorporating a twist-averaging 
procedure in our simulations. For the lowest tempera¬ 
ture, 0 = 0.0625 and = 1, FSC(a) leads to reasonable 
agreement with analytical approximations (see below), 
and smoothly connects to the RPIMC data, for Vg ^ 5, 
cf. Fig. For smaller Vg and higher 0 the formula is 
not applicable. On the other hand, the FSC of Brown 
produces energies that are systematically too high and 
has to be discarded [53]. Due to the lack of applica¬ 
ble high-density FSC, we performed additional CPIMC 
simulations for a broad range of particle numbers up to 
Nmax = 400, allowing for an extrapolation to the macro¬ 
scopic limit, for Vg = 0.1 and 0 = 0.0625; 0 = 0.5. For 
0 = 0.0625 an additional point at Vg = 0.3 could be ob¬ 
tained [53], cf. the diamonds in Fig. These accurate 
data will be a suitable starting point for the construction 
of FSC formulas that are applicable at high densities. 

To obtain independent analytical results for the macro¬ 
scopic UEG, we now compute the exchange correlation 
energy including, in addition to Hartree-Fock [32], the 
two second order diagrams (Montroll Ward and e"*) [50] . 
The two results [cf. Fig. [^, converge for low r^, even¬ 
tually reaching the Hartree-Fock asymptote (horizontal 
line). For rg > 0.1, MW and start to deviate from one 
another, and we expect the exact result to be enclosed 
between the two m- Reliable predictions are possible up 
to rg - 0.8, for 0 = 0.0625, and r, - 0.55, for 0 = 0.5 



FIG. 5: (Color online) Exchange correlation energy (times Vg) 
of the macroscopic polarized UEG at 0 = 0.0625 (blue) and 
0 = 0.5 (red). Open symbols: RPIMC results [27]. CPIMC 
(a): our results with FSC from Ref. [52] ■ RCPIMC (b): 
our data (3 points) with numerical extrapolation N ^ oo 
[54]. Dashes: analytical e'*-approximation [50], dots: fit of 
Ref. [32] ■ For better visibility the curves for 0 = 0.5 are 
upshifted by 0.2. 

[50] . In Figurewe also include the fit of Ref. [32] that 
shows, overall, a very good behavior, but is too low at 
Tg —>■ 0, with the deviations growing with 0 m- 

To summarize, we have presented first-principle config¬ 
uration PIMC results for the UEG at finite temperature 
that have no sign problem at high degeneracy, rg < 0.4, 
and allow for reliable predictions up to = 1. This 
makes CPIMC a perfect complementary approach to di¬ 
rect fermionic PIMC and to RPIMC that cannot access 
high densities, and our results indicate that the previ¬ 
ous RPIMC data are not reliable in the range rg < 4. 
The present results will be important for dense quan¬ 
tum plasmas at finite temperatures, that are relevant for 
warm dense matter, in general, and for ICF, in particu¬ 
lar. Since here the electrons are typically unpolarized we 
tested our CPIMC approach for this case. Although the 
sign problem is more severe than for the polarized situa¬ 
tion CPIMC is well capable to produce very accurate ab 
initio hnite temperature results that smoothly connect 
to available T = 0 data [50]. A systematic analysis is 
presently under way. 

Finally, the obtained accurate results for the exchange- 
correlation energy provide benchmarks for finite temper¬ 
ature DFT, RPIMC [2Z], novel independent QMC sim¬ 
ulations [311EZ] and analytical fits [32|. Although the 
fermion sign problem is not removed by our approach, 
the proposed combination of CPIMC with DPIMC (or 
RPIMC) provides, for the UEG, a practical way around 
this problem. 

We acknowledge stimulating discussions with V. Fili- 
nov, J.W. Dufty, V.V. Karasiev and S. Trickey and E. 
Brown for providing information on the FSC (a) used 
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